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Abstract. Electrostatic actuators are simple but important switching devices 
for MEMS applications. Due to the difficulties associated with the electrostatic 
nonlinearity, precise mathematical description is often hard to obtain for the 
dynamics of these actuators. Here we present two sharp theorems concerning 
the dynamics of an undamped electrostatic actuator with one-degree of free- 
dom, subject to linear and nonlinear elastic forces, respectively. We prove that 
both situations are characterized by the onset of one-stagnation-point periodic 
response below a well-defined pull-in voltage and a finite-time touch-down or 
collapse of the actuator above this pull-in voltage. In the linear-force situation, 
the stagnation level, pull-in voltage, and pull-in coordinate of the movable elec- 
trode may all be determined explicitly, following the recent work of Leus and 
Elata based on numerics. Furthermore, in the nonlinear-force situation, the 
stagnation level, pull-in voltage, and pull-in coordinate may be described com- 
pletely in terms of the electrostatic and mechanical parameters of the model 
so that they approach those in the linear-force situation monotonically in the 
zero nonlinear-force limit. 



1. Introduction 

It is well known that electrostatic actuators play a basic and important role in 
microelectromechanical systems (MEMS) technology as one of the most commonly 
used microcircuit components and mathematical modeling and numerical simula- 
tion are essential for understanding various qualitative and quantitative properties 
of such devices [SJ[TT]. However, due to the inverse-square type nonlinearity in the 
governing equations associated with the Coulomb electrostatic forces, it is a chal- 
lenging task for analysts to gather a precise picture of the mathematical behavior of 
the systems and numerical simulations have been resorted to in the investigations 
of the modeling of actuators. So far, overwhelmingly, mathematical studies have 
been concentrated on equations governing steady state solutions [21 HI M 033 H5 
or their parabolic extensions [H [T^] and little work has been done concerning the 
dynamical behavior of the problem which is of a wave-equation nature and of ob- 
vious importance. (In [TJ [BJ, some existence, uniqueness, and continuous depen- 
dence results are obtained for the evolution problems related to the MEMS wave 
equations. We thank the referee for informing us these developments.) In this pa- 
per, as a first attempt, we shall study the dynamics of an undamped actuator of 
one-degree of freedom described by a second-order ordinary differential equation, 
as a simplified model arising from the original spatially dependent wave equation 
description [IT]. The simplicity of the model enables us to obtain two sharp the- 
orems concerning the dynamics of the actuator, subject to linear and nonlinear 
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elastic forces, respectively, which roughly say that in both situations the dynamic 
response of the actuator under an applied voltage V is characterized by the onset of 
one-stagnation-point periodic motion when V is below a well-defined pull-in voltage 
and a finite-time touch-down or collapse of the actuator happens when V is above 
this pull-in voltage. More precisely, in the linear-force situation, the stagnation 
level, pull-in voltage, and pull-in coordinate of the movable electrode may all be 
determined explicitly, following the recent work of Leus and Elata [8] based on nu- 
merics, while, in the nonlinear-force situation, the stagnation level, pull-in voltage, 
and pull-in coordinate may be described completely in terms of the electrostatic 
and mechanical parameters of the model so that these quantities approach those 
in the linear- force situation monotonically in the zero nonlinear- force limit. It is 
interesting that, under a fairly natural convexity assumption, our method for the 
linear-force situation is completely applicable to the nonlinear situation so that an 
understanding of the global dynamics of the undamped actuator is achieved. 

The rest of the paper is organized as follows. In Section 2, we review the physical 
description and equation of motion of the actuator subject to a linear elastic force 
and the work of Leus and Elata [8], which motivates our study to follow. In Section 
3, we use the energy conservation law as a first integral to establish the existence of 
a stagnation time and stagnation level for the motion when the applied voltage is 
below a critical value. In Section 4, we use the results of Section 3 to show that the 
motion with a stagnation point actually gives rise to a symmetric periodic motion 
of the actuator. In Section 5, we show that when the applied voltage is above the 
critical value identified in [8] and Section 4, the solution is non-periodic and leads 
to a touch-down or collapse situation of the device. In Section 6, we summarize our 
study as a theorem which describes the full dynamical behavior of the actuator and 
confirms the results stated in [8], when the elastic force is linear. In Section 7, we 
extend our method to study a more general model when the elastic force is cubicly 
nonlinear which is widely useful in applications (cf. [5] and references therein). 
We conclude by indicating how our method may be used to tackle models under 
more general nonlinear elastic forces. In Section 8, we spell out some simple but 
generic conditions under which the dynamic touch-down or collapse phenomenon 
occurs universally as a result of a sufficiently high applied voltage. 



2. Physical model and governing equation 

Following [5], the parallel-plate actuator of one-degree of freedom under study 
consists of a top electrode of mass m and area A that is suspended on a linear elastic 
spring with the Hooke constant k, above a fixed bottom electrode that is coated 
with a dielectric layer of thickness do > 0. The initial gap or distance between the 
top electrode and the dielectric is g > 0. The fixed bottom electrode is electrically 
grounded and a step-function type voltage will be applied to the top electrode so 
that it is zero before time t = and assumes the constant value V > when t > 0. 
Such an application results in the one-dimensional motion of the top electrode away 
from its initial equilibrium position with the coordinate x = to that with x > 
so that the Hamiltonian of the motion, which is the sum of the kinetic energy, 
elastic potential arising from a linear elastic force, and electrostatic potential of the 
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deformable capacitor subject to the voltage V, is given by the expression [8] 

(2.1) H = -mx 2 + -kx 2 - . £oAV2 . . 

2 2 o / d <> 



where x = dx/dt is the velocity of the moving top plate, eo the permittivity of the 
free space, and £r the relative permittivity of the dielectric. 

In order to simply the notation, we may conveniently adopt the normalized 
variables and parameters |8, as follows 

/■„ ~s x H „ jk sr,AV 2 Tr9 dn 

(2.2) -^x, jr^^H, J-t^t, -JL^^y 2 , -2- = £. 

g kg^ V m kg a ge^ 

Hence the Hamiltonian (|2.ip takes the form 

, s 1 9 1 9 1 V 2 

^ 2 + 2 ~ 2 (iTT 3 ^) ' 

As a consequence, the momentum equation governing the motion of the top plate- 
shaped electrode is [8] 

9H IV 2 

(2 - 4) ^~^-^ + 2(eTi-^' ^ >0 - 

Since the actuator is free before the onset of the positive voltage V > starting 
at t = 0, the initial condition for the coordinate variable of the top electrode is 
given by 

(2.5) x(0) = 0, i(0) = 0. 

The simple initial value problem consisting of (|2.4I) and (12 .5[) and describing the 
dynamic response of an undamped actuator [8J, subject to a step- function type 
applied voltage and a linear elastic force, will be our main focus in the first part of 
our study here. Although the nonlinearity arising from the Coulomb electrostatic 
force appears to be singular at x = l + the touch-down or collapse actually occurs 
at 

(2.6) x = l, 

due to the fact that the bottom electrode is coated with a dielectric layer of thickness 
do > and relative permittivity e r , which results in an effective insulator layer of 
the 'thickness' £ > given in (|2.2p in the normalized units. 

Based on numerical simulations, Leus and Elata [S] find that there is a critical 
voltage, called the 'dynamic pull-in voltage', Vj^pj, so that below Vjjpj the dynamic 
response of the actuator is periodic, above Vfjpj, the response is non-periodic so 
that the motion is characterized by a positive velocity which results in a finite-time 
touch-down or collapse, and at the critical voltage Vj^^ the movable electrode 
converges to an unstable equilibrium position as t — > oo, called the 'dynamic pull- 
in' position, ^Tjpj > 0. In order to calculate the dynamic pull-in voltage and 
dynamic pull-in position, Leus and Elata [8] investigated the periodic solutions 
which, according to their numerical simulations, would oscillate between the initial 
state x = and the state of the maximum displacement, x$, where the velocity 
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vanishes, called the stagnation displacement. In order to find xs, they use the 
energy conservation law 

(2.7) fr| t>0 = fr| t=0 = _I^_ 

to get the first integral 

(2.8) i 2 +x 2 -- - = 

1 ' (£ + 1 -*)(£ + 1) 

for any solution at t > 0. Hence, the stagnation position Xs may be obtained by 
setting x = in (12.81) to yield 

(2.9) (£ + l)xl - (£ + l) 2 ^ + V 2 = 0. 

These equation has two positive roots, one is below (£ + l)/2 and the other above 
(£ + l)/2; both are below £ + 1. Leus and Elata dismiss the larger root as non- 
physical and choose the smaller root to determine the stagnation position, 

r9in x _ (e + i)W(g + i) 4 -4(g + i)^ 2 " 

xs- 2(e + 1) 

The formula (12.10)) clearly indicates that the stagnation position increases with 
regard to the applied voltage until it fails to exist when the argument in the square 
root of (|2.10p becomes negative. As a consequence, they conclude that the dynamic 
pull-in voltage and dynamic pull-in position are given by the explicit formulas 

(2-11) *t)pi=~(l + ! , *Dpi=~(l+0- 

In the first part of this paper consisting of Sections 3-6, we aim at establishing 
the above global dynamic picture for an undamped parallel-plate actuator subject 
to a linear elastic force rigorously. In the second part of the paper (Section 7), we 
will use the method developed in detail in the previous sections to obtain the same 
dynamic picture when the actuator is under a cubicly nonlinear elastic force. In 
doing so, we will see that our method may readily be extended to investigate more 
general models. 

3. Use of the first integral 
To proceed, we rewrite the conservation law (12.81) as 

We see from (I3.1[) that any solution x(t) of the initial value problem consisting of 
(|2.4p and (|2.5[) must remain nonnegative-valued for all t, 

(3.2) x(t) > 0. 

For convenience, we denote the right-hand side of (13. ip by f(x). Then it takes 
the factored form 

«"» '« - «tM- - « + * + (i^i)) - wrhrf^ 

We see that, when the voltage V is low such that 
(3.4) y2 < l (e + 1)3j 
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the quadratic polynomial q(x) has two real roots, say xx and x 2 , given by the 
expressions 



1 1 / W 2 

(3-5) *i,2 = ^ + l)T^ (£ + l) 2 - " 



2^ ' ' 2^^ ' (£+1) 
It is important to observe the property 

(3.6) < x x < x 2 < £+ 1, 
which allows us to express fix) as 

(3.7) f(x) = j—^—rixt- x)(x 2 - x). 

As a consequence, we see that the only physically acceptable solution has to satisfy 
the boundedness condition 

(3.8) < x(t) < xx, 

otherwise it would violate the non- negativity of the right-hand side of (|3.ip . 
With (I3.8[) in mind, we decompose (|3.ip into two pieces, 



(3.9) x = J ^ - - _ ^ (xx - x)(x 2 - x), t>0, 

and 



(3.10) x = - J ^ - | - ^ (xx ~ x)(x 2 ~ x), t>0, 

which may be viewed as two first integrals of (|2.4[) . It is a comfort to note that 
the initial condition x(0) =0 automatically implies the desired initial condition for 
velocity, x(0) = 0, stated in (|2.5[) . 

Using the condition x(t) > 0, we see that the initial condition x(0) = selects 
the equation (|3.9p over (|3. 10[) for t near t — 0. In other words, a solution always 
starts its evolution following (|3.9[) rather than p. 101) . 

Let x(t) be a local solution of (|3.9p starting from x(0) = 0. It will be useful to 
estimate time t = tx > needed for a;(t) to cross the level x = xx/2 (say). To this 
end, we use 

(3.11) 0<a;(t)<y, 0<t<tx 



and (|3.9p to arrive at the inequality 



(3.i2) ^^wttA^-t)' °^^ h - 

Integrating ()3.12p over < t < tx, we find the following upper bound for tx, 



(3.13) tx<2V2~J J + 1 . 

V ZX2 — Xx 

We next prove that there is a later time ts > tx such that x(t$) = Xx- In [8], 
ts is called the stagnation time corresponding to the stagnation position xs = x\ 
since the velocity of the motion vanishes at t = ts when the top electrode stretches 
to the position x = x$. 
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Indeed, for t > t\ so that x(t) stays below xi, that is, 
(3.14) ^ < x(t) < x x (= x s ), 



we obtain from (|3.9p the differential inequality 

(3-15) i > y 2 ^ ^y(^2 - ~ x, t > t x . 

Integrating (|3.15l) . we have 



(3.16) 2\J^-y/ Xl - x(t)j > ^ (x 2 - xi)(t - h), t>ti. 

Such a relation establishes the expected property that there must be a finite ts > t\ 
such that 

(3.17) x(t s ) = x s = xx 
as claimed. 

4. Periodic solutions 

We consider what happens after the stagnation time ts- Naturally, we may intend 
to use (|3.9[) again to extend the solution. However, this cannot be done because 
(I3.9[) results in an increasing solution but the right-hand side of (I3.9[) is invalid for 
x > Xs = X\. Therefore, we need the solution to decrease after < s - Fortunately, 
such a feature is reflected in (I3.10[) . In other words, we are led to use (|3.10l) to 
extend the solution after ts- In fact, we can do so explicitly by setting 

(4.1) x(t) -> x(2t s - t), t>t s . 

In other words, the definition 

x(t), 0<t<t s , 
x(2t s -t), t s <t<2t s 

gives us a solution of (I2.4[) which satisfies (|3.9[) for < t < ts and (|3.10l) for 
ts < t < 2ts- This construction clearly produces a periodic solution with period 

(4.3) t p = 2t s . 

In view of <|3.13p and (|3.16|1 . we arrive at an upper bound estimate for ts, 



(4.2) x(t) 



(4.4) ia < 2(( i±i^i + ^/ « + > 



X2 — X\ V 2^2 — X\ 

5. Non-periodic or touch-down solution 



We now study what happens when (J3T4J) is violated, 
(i) We first assume that 

(5.1) y 2 >i(£ + l) 3 . 

Since (|5.1[) implies that the quadratic function q{x) defined in Section 3 is positive 
definite, 

(5.2) mm{q(x)} = " ^ a 2 > 0. 
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-q(x) > 



Hence (|3.9[) leads to 

(5 - 3) 1 ve+i 

Integrating (|5.3p over [0,t], we find 
(5.4) 2v^(t)> 



fx, t > 0. 



which indicates that there must be a finite time ic > 0, called the "contact time", 
such that 



(5.5) 



x(tc) — lim x(t) = 1. 



Beyond tc, the solution does not exist and (|5 . 5[) since now the top electrode is 
in contact with the dielectric layer. In particular, when £ = so that the bottom 
electrode is not coated with dielectric, (15.51) implies that the switch device is "short- 
circuited" . Thus, the critical level 



(5.6) 



V 2 



1 



-(£ + l) 3 



which borders the dynamic touch-down and periodic response situations of the 
device, is rightfully referred to as the pull-in voltage [5]. 
(ii) We now assume (|5.6p . From (|3.5[) . we see that 



(5.7) 

and (|3.9p becomes 
(5.8) 



q(x) 



itt- 



1) 



Since £ > 0, we see that £ + 1 > (£ + l)/2 and the right-hand side of (|5.8[) is 
Lipschitzian around x = (£+l)/2. As a consequence, the uniqueness of the solution 
to the initial value problem of the ordinary differential equation (|5.8I) holds in a 
neighborhood of x = (£ + l)/2. Since ir = (£ + l)/2 is an equilibrium point of 
(|5.8p . it is not attainable at a finite time i > by a solution different from this 
equilibrium point. Consequently, we must have 



(5.9) 



lim x{t) = i(£ + l), 



which is the limiting position coordinate of the top electrode at the pull-in volt- 
age given in (I5.6[) and referred to as the pull-in position coordinate [5] of the top 
electrode. 



6. A PRECISE DYNAMIC RESPONSE THEOREM 

In summary of the study carried out in the previous sections, following [S], a 
complete description of the dynamic response of an undamped electrostatic MEMS 
actuator subject to a linear elastic force is obtained as follows. 

Theorem 6.1. Consider the differential equation {2.1$ governing the dynamics 
of an electrostatic actuator with a top movable electrode subject to a linear elastic 
force and the free initial condition &2.5\) . When the applied voltage is low such that 
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V < i(£ + l) 3 , the response of the top electrode is periodic characterized by a 
stagnation time t$ > at which the electrode stretches to its maximum distance 

and the electrode follows the same route to return to its original free initial position 
att = 2t s . In other words, the response is represented by a period tp — 2t s function 
which is symmetric about the point t — t s . When the applied voltage is high such 
that V > 1) 3 / 2 , the response of the top electrode is such that, it is driven into 

contact with the dielectric layer coated over the bottom electrode plate or the bottom 
electrode plate itself if it is not coated, at a finite time t c > 0. In other words, 
x{t) — > 1 as t — > t c . As a consequence, the critical pull-in voltage is identified 
explicitly as Vjjp^ — |(£ + l) 3 / 2 . When the applied voltage is at the level Vjj^, 
the plate-shaped electrodes will never make contact with each other but the top 
electrode approaches an equilibrium position, or the pull-in position, characterized 
by x(t) — > xjjpj = |(£ + 1) as t — > 00. 

In the next section, we will see how the method developed in proving Theorem 
16.11 for a linear-elastic force situation can be adapted to obtain an equally sharp 
dynamic response theorem for a widely useful nonlinear-elastic force situation. 

7. A NONLINEAR ELASTIC FORCE SITUATION 

We now consider the undamped parallel-plate actuator subject to a nonlinear 
elastic force. A typical model is given by the Hamiltonian ^ [5] 

(7.1) H = \mx 2 + \kx 2 + h 3 x 4 S ° AV2 



(in 



where k 3 > is usually a small constant which gives rise to a cubicly nonlinear 
clastic force. Using the new variables and parameters following (|2.2j) and setting 

(7.2) K= __ 

the Hamiltonian (|7.1[) is normalized to take the form 

(7.3 H = -x 2 + -x 2 + -kx a -, 

y ' 2 2 4 2 (£ + 1 - x) ' 

so that the momentum equation governing the motion of the top electrode reads 

OH 3 1 V 2 

~ — — X KX ~\ — ~, , 

dx 2 (£ + 1 - xY 

As before, in view of the energy conservation, we can obtain the first integral of 
CLU) as 

(7.5) ± 2 = ^ (J^-(£ + i- x \ x -1(£ + 1 - x ) x z 

For the normalized coated dielectric layer of thickness £ > 0, the solution of interest 
satisfies 

(7.6) < x < 1. 



(7.4) x = — — = -x - kx 3 + - — ■ t > 0. 
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Thus, wc need to require 

(7.7) g(x) = -(£ + i- x ) x -!l(£ + i- X )x 3 > 

for consistency in (|7.5I) and (|7.6I) . Motivated by the study in the previous sections 
where k = 0, we shall maintain the global strict convexity of g(x) by imposing 

(7 - 8) K < 3T?Ti?' 

which is compatible with the fact that k% is small (cf. (I7.2[l ) in applications. 

We need to investigate when the global minimum of g{x) becomes negative. We 
can check that 

(7.9) ff'(0) = -(C + l)<0 > ff'(£ + l) = (e + l)(f(£ + l) 2 + l) >0. 
Hence there is a unique point xq £ K such that 

(7.10) </(a;o) = 2kxI - ~k{£ + l)x 2 + 2x - (£ + 1) = 0, < z < £ + 1, 

since g(x) is strictly convex due to (I7.8[) . Such an xo is the global minimum point 
of the function g(x). It is interesting to note that xq depends only on £ and k but 
not on V, which enables us to treat the following three cases individually in a nice 
way. 

Case (i): V is sufficiently low so that g(xo) < 0. Since it is clear that 

(7.11) g(x)>0, x e (-oo,0] U [£ + l,oo), 

we see that there are exactly two points, say x\ and x^, such that 

(7.12) g(x* 1 )=g(x* 2 ) = 0, x\ e (0,^), e (x ^ + 1). 
Therefore g(x) enjoys the factorization 

(7.13) g(x) = (xl -x)(xl -x)q(x), Vx, 

where g(a;) is a positive- valued quadratic polynomial. Inserting 13[) into (|7.5I) . 
we arrive at 

(7.14) i2 = TTl (xt - x){x* - x)q{x). 

Consequently, we may duplicate the methods in Sections 3 and 4 to show that there 
is a finite time t$ such that a solution satisfying the initial condition (I2.5[) climbs 
to the stagnation level x§ = x\ when t = tg, 

(7.15) ^)=x*=^, 

and descends back to the initial position x — after the same amount of elapse 
of time, tg. In other words, we obtain a symmetric periodic solution of period 

It will be interesting to compare x§ with xs obtained in Section 4, when V < 
Vrjpj. In fact, with the notation of Section 4, we have 

(7.16) g{ Xl ) = -^ + l- Xl )xl<0, 



(7.17) g'{x x ) = -^ + l)-(l + KX l)J^ + l)2-^- <0 . 
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From (|7.17|) and g"{x) > 0, we have x\ < xo; from (|7.16|) and g'(x) < for 
x S [0, xo), we conclude that 



(7.18) xS=x* 1 <x 1 =x a = ±(t + l)-± ) J(£+ l) 2 - j^- y 

whenever V < Vq^. That is, in the situation of a nonlinear elasticity for the top 
electrode, the stagnation distance is shorter than that of the linear situation. 

Case (ii): V is sufficiently high so that g(xo) > 0. We may use the method in 
Section 5 to show that there is a finite contact time, such that 

(7.19) x{t* c ) = lim x{t) = 1. 

In other words, the top electrode plate stops at time t = i£ to take a touch-down 
contact position with the bottom electrode plate. 

Case (iii): V assumes a critical value, say Vc, so that g(xo) = 0. That is, xo 
satisfies 

(7.20) = (£ + l-x )x + |(£ + l-x )xg. 

Since Xo also satisfies (I7.10[) . it is a double root of the equation g = 0, which implies 
the factorization 

(7.21) g(x) = {x -x) 2 q{x), 

where q(x) is a positive- valued quadratic polynomial. Substituting (I7.21[) into (|7.5p . 
we have 



(7.22) i= L^)_\ X0 _ xl 

y £ + 1 - x 

Again, we may use the method in Section 5 to show that a solution of (|7.22|) 
satisfying the initial condition (|2.5j) is defined for all t > and has the monotone 
limit 

(7.23) lim x(t) = x . 

t— >QO 

Following Section 5, we define the critical voltage Vq given in (|7.20[) to be the 
dynamic pull-in voltage, ^Qpp that is 



(7.24) = ^ + 1) ((£ + 1 - xo)x + f (£ + 1 - ^o)x3) , 

and the corresponding asymptotic position coordinate xq, given in (|7.10p . to be the 
dynamic pull-in position coordinate, xj^pj, 

(7.25) xj5 pi = Xq. 
However, since 

(7.26) g'(x m ) = g' Q(l + £)) = K + I) 3 < 0, 

and g"(x) > 0, we see that there holds 

(7.27) x = x^ pi > x Dp j, if k > 0. 
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Furthermore, use h(x) to denote the function 
(7.28) h(x) = (t + l-x)x + ^ + l-x)x\ 

so that the right-hand side of (|7.20l) is simply h(xo). Then h'(x) — —g'(x). In 
particular, h'{x) > for x < Xq, which implies that h(xjy ^) > M^Dpi) - ^ n 
immediate consequence of this fact, combined with (|7.20p . is that 

„ 29) <g£>!|* + « ({ + 1) . lt „ >0 , 

in particular, 

(7.30) ^Dpi>^D P i> if*>0, 

which is a naturally expected result, which says that the presence of a cubicly 
nonlinear elastic force for the motion of the movable electrode of the actuator 
enhances both dynamic pull-in voltage and dynamic pull-in position coordinate. 

It is worth noting that, in the limiting situation when n = 0, the formulas (|7.10p 
and (|7.24[) recover the linear elasticity results obtained earlier, so that 

(7.31) ^.^^ = !(£ + !), ^ pi ^F Dpi = i(£ + l)i, k^O. 

The study just carried out enables us to arrive at the following second sharp 
dynamic response theorem concerning the MEMS actuator subject to a cubicly 
nonlinear elastic force. 



Theorem 7.1. Consider the differential equation J7.^| ) describing the dynamics 
of an undamped electrostatic actuator subject to a cubicly nonlinear elastic force 
and the free initial condition \2. 5\) . governed by the normalized Hamiltonian Ij7.3\ ), 
and assume that the nonlinear elasticity is weak so that J7.<$[ ) holds. Then the 
full picture of the dynamic response of the system is completely determined by the 
unique solution xj^ ^ — Xq of the cubic equation which necessarily lies in 

the interval (0,£ + 1) and defines the dynamic pull-in voltage Vjj ^ through the 

formula \1.2J$ . When the applied voltage is low such that V < Vjj p the response 

of the actuator is periodic characterized by a stagnation time t* s > at which the 
top movable electrode stretches to its maximum distance x* s so that 

1 1 / W 2 

(7.32) 4 <a;s = _(£ + i)__J(£ + i)2___ if V <V Dpi , 

and the top electrode follows the same route to return to its initial position at 
t = 2t* s . In other words, the response is represented by a period tp = 2t* s function 
which is symmetric about the point t = t* s . Furthermore, as a function of k and 
V , the stagnation position coordinate x* s is smooth, strictly decreasing with respect 
to k, and increasing with respect to V . When the applied voltage is high such that 
V > Vjj the response of the actuator is so that, the applied voltage drives the top 

electrode plate into contact with the dielectric layer coated on the bottom electrode 
plate or the bottom electrode plate when there is no coated dielectric layer, at a 
finite time t c > 0. In other words, x(t) — > 1 as t — > t c . When the applied voltage is 
at the level Vj^ ^, the top electrode plate will never make contact with the dielectric 
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layer or the bottom electrode plate but approach an equilibrium position, or the pull- 
in position, characterized by x(t) — > xjj^ as t — > oo. Moreover, x*j^ ^ and Vjjp^ 
depend on n smoothly and are strictly increasing, so that they converge to Xj^^ and 
Vjjpi in the zero nonlinear- elastic force limit k — > 0. 

Proof. We only need to prove the smoothness and monotonicity of Xg , Xj~^^ , and 
Vgp- with respect to the designated variables. 

To proceed, assume V < Vjj-pi ^ n ^ n * s situation, we first notice that the condi- 
tion (|7.8|) ensures that g'(xg) = g'(x*) < since g'(x ) = and g"{x) > for all x. 
Thus, the implicit equation 

(7.33) G{x\, k, V) = ^ - (C + 1 - xl)x\ - | (£ + 1 - x\){xlf = g{xl) = 
gives us 

dG 

(7.34) — (x* 1 , K ,V)=g'(x* 1 ) < 0, 

and the smooth dependence of x\ on k and V follows as a result of the implicit 
function theorem. Besides, after implicit differentiation, we have 

dG\dx* 1 l ft . .,,.. 3 



(7-35) {—J^ = -(*+!-*;)(*;)' >0, 

/ dG \ dx\ 2V 

^ [WJW = -?TT <0 - 

Therefore, combining (|7.34j) - (|7.36j) . we arrive at the conclusion 

(7.37) M <0 , M >0 . 

We now show that xjjpj an d ^rjpj are smoothly defined increasing functions of 
the parameter k. To this end, recall that a^p- = xq satisfies the equation (|7.10|) . 
which may be rewritten as F(xo, k) — 0. Of course, 

(7.38) £=^°-^ + 1 ^ + ^) >0 ' 

due to the condition (17.81) . which establishes that xq is a smooth function of k 
determined implicitly by F(xq, k) = 0. Using (|7.10[) . we have 

8F\ dx fdF\ f n ^ 3.. a 



axo / d« \ ok / \ 2 

(7.39) = I (^0 - + 1)J >0, 

in view of (|7.27l) and the fact that xpjpj = + 1) obtained earlier. Combining 
([?T38l) and (|7T39| . we have > as claimed. 

With the notation ([7^5]) . we rewrite ([7^4]) as Vq ■ = + l)h(x ). Therefore 

(7-40) = '[^±4 ( h'(x )^ + k + 1 - > 0, 



d« 2 V h(xo) \ dn 2 

since h'(xo) — —g'(xo) = and xq < £ + 1 
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The proof of the theorem is complete. □ 

The physical meaning of the monotone dependence of Xg = x\ on V and K is 
natural to expect since a larger voltage V drives the top electrode plate down further 
but stronger elasticity characterized by a larger k parameter gives greater stiffness, 
which in turn leads to a smaller stretching distance for the periodic motion. 

We also note that our method may be used to study the model governed by the 
following normalized extended Hamiltonian 

1 IV 2 
(7.41) H = -x 2 + $(x) 



2 w 2(£+l - x y 

where $(cc) > with $(0) = is a general elastic potential density function and 
damping is absent. In fact, since the energy conservation law gives us the relation 

(7 - 42) l2 = FT)FM- 2$(l) ' 

we see that the crucial condition to be imposed so that our method works would 
be the strict convexity of the function 

V 2 $(x) 
(7.43) *( x ) = *-2(£ + l-x)-fl-, x>0. 

8. Touch-down as a universal dynamical property 

In this section, we shall identify some simple conditions under which the touch- 
down phenomenon inevitably takes place when a coupling parameter resembling the 
applied voltage is sufficiently large. For greater generality, we will include damping 
in the model so that the normalized equation of motion assumes the form 

(8.1) x + fix + f(x,t) = Xg(x,t), 

where fi > is a damping constant, a physically meaningful solution x = x(t) stays 
in the interval [0,a] for some a > 0, and / and g are functions of x and t so that 
there are positive constants C\ and C2 such that the bounds 

(8.2) sup{|/(x,i)||xe [0,a],t>0} < d, 

(8.3) mi{g(x,t)\x£ [0,a),t>0} > C 2 . 

hold. 

We will assume a to be the touch-down position for the solution so that the 
function g, a term resembling the Coulomb force, may be singular at x = a. We 
have 

Theorem 8.1. If X > is sufficiently large, there is an associated finite time 
t c > 0, depending on X and a, such that the solution of i8.1\) subject to the initial 
condition h2. 5\) monotonically climbs to the touch-down position at t — t C - That is, 

(8.4) x(t c ) = lim x(t) = a and x(t) > Vt£(0,t c ). 

t-¥tc 

Proof. First assume fi > 0. In view of (|8.1| - (|8.3|) and the initial condition (|2.5|) . 
we have, after multiplying (|8.1j) by e M * and integrating and using the condition 
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< x(t) < a, the relation 




(8.5) 



> -(AC 2 -d)(e^-l), 



which establishes x[t) > for all i > provided that AC2 > C\ and leads to the 
lower estimate 



Hence there is a finite time tc > such that (|8.4|) holds when AC2 > C\ . 

Next, if n = 0, then an integration of (|8.1|) simply gives us x(t) > (AC2 — C\)t 
as before and the rest follows obviously. 



In summary, we have presented two theorems which give precise descriptions 
of the global dynamic behavior of an electrostatic undamped actuator, subject 
to linear and cubicly nonlinear elastic forces, arising in MEMS technology. The 
method developed is applicable to the study of the models with more general elastic 
forces. 
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(8.6) 




Thus the theorem is proven. 



□ 
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